Radiomic Phenotypes for Improving Early Prediction of Survival in Stage III Non-Small Cell Lung Cancer Adenocarcinoma after Chemoradiation.

We evaluate radiomic phenotypes derived from CT scans as early predictors of overall survival (OS) after chemoradiation in stage III primary lung adenocarcinoma. We retrospectively analyzed 110 thoracic CT scans acquired between April 2012-October 2018. Patients received a median radiation dose of 66.6 Gy at 1.8 Gy/fraction delivered with proton (55.5%) and photon (44.5%) beam treatment, as well as concurrent chemotherapy (89%) with carboplatin-based (55.5%) and cisplatin-based (36.4%) doublets. A total of 56 death events were recorded. Using manual tumor segmentations, 107 radiomic features were extracted. Feature harmonization using ComBat was performed to mitigate image heterogeneity due to the presence or lack of intravenous contrast material and variability in CT scanner vendors. A binary radiomic phenotype to predict OS was derived through the unsupervised hierarchical clustering of the first principal components explaining 85% of the variance of the radiomic features. C-scores and likelihood ratio tests (LRT) were used to compare the performance of a baseline Cox model based on ECOG status and age, with a model integrating the radiomic phenotype with such clinical predictors. The model integrating the radiomic phenotype (C-score = 0.69, 95% CI = (0.62, 0.77)) significantly improved (p<0.005) upon the baseline model (C-score = 0.65, CI = (0.57, 0.73)). Our results suggest that harmonized radiomic phenotypes can significantly improve OS prediction in stage III NSCLC after chemoradiation.

Leveraging machine learning predictive biomarkers to augment the statistical power of clinical trials with baseline magnetic resonance imaging.

A key factor in designing randomized clinical trials is the sample size required to achieve a particular level of power to detect the benefit of a treatment. Sample size calculations depend upon the expected benefits of a treatment (effect size), the accuracy of measurement of the primary outcome, and the level of power specified by the investigators. In this study, we show that radiomic models, which leverage complex brain MRI patterns and machine learning, can be utilized in clinical trials with protocols that incorporate baseline MR imaging to significantly increase statistical power to detect treatment effects. Akin to the historical control paradigm, we propose to utilize a radiomic prediction model to generate a pseudo-control sample for each individual in the trial of interest. Because the variability of expected outcome across patients can mask our ability to detect treatment effects, we can increase the power to detect a treatment effect in a clinical trial by reducing that variability through using radiomic predictors as surrogates. We illustrate this method with simulations based on data from two cohorts in different neurologic diseases, Alzheimer's disease and glioblastoma multiforme. We present sample size requirements across a range of effect sizes using conventional analysis and models that include a radiomic predictor. For our Alzheimer's disease cohort, at an effect size of 0.35, total sample size requirements for 80% power declined from 246 to 212 for the endpoint cognitive decline. For our glioblastoma multiforme cohort, at an effect size of 1.65 with the endpoint survival time, total sample size requirements declined from 128 to 74. This methodology can decrease the required sample sizes by as much as 50%, depending on the strength of the radiomic predictor. The power of this method grows with increased accuracy of radiomic prediction, and furthermore, this method is most helpful when treatment effect sizes are small. Neuroimaging biomarkers are a powerful and increasingly common suite of tools that are, in many cases, highly predictive of disease outcomes. Here, we explore the possibility of using MRI-based radiomic biomarkers for the purpose of improving statistical power in clinical trials in the contexts of brain cancer and prodromal Alzheimer's disease. These methods can be applied to a broad range of neurologic diseases using a broad range of predictors of outcome to make clinical trials more efficient.

Interactive Machine Learning-Based Multi-Label Segmentation of Solid Tumors and Organs.

We seek the development and evaluation of a fast, accurate, and consistent method for general-purpose segmentation, based on interactive machine learning (IML). To validate our method, we identified retrospective cohorts of 20 brain, 50 breast, and 50 lung cancer patients, as well as 20 spleen scans, with corresponding ground truth annotations. Utilizing very brief user training annotations and the adaptive geodesic distance transform, an ensemble of SVMs is trained, providing a patient-specific model applied to the whole image. Two experts segmented each cohort twice with our method and twice manually. The IML method was faster than manual annotation by 53.1% on average. We found significant (p < 0.001) overlap difference for spleen (DiceIML/DiceManual = 0.91/0.87), breast tumors (DiceIML/DiceManual = 0.84/0.82), and lung nodules (DiceIML/DiceManual = 0.78/0.83). For intra-rater consistency, a significant (p = 0.003) difference was found for spleen (DiceIML/DiceManual = 0.91/0.89). For inter-rater consistency, significant (p < 0.045) differences were found for spleen (DiceIML/DiceManual = 0.91/0.87), breast (DiceIML/DiceManual = 0.86/0.81), lung (DiceIML/DiceManual = 0.85/0.89), the non-enhancing (DiceIML/DiceManual = 0.79/0.67) and the enhancing (DiceIML/DiceManual = 0.79/0.84) brain tumor sub-regions, which, in aggregation, favored our method. Quantitative evaluation for speed, spatial overlap, and consistency, reveals the benefits of our proposed method when compared with manual annotation, for several clinically relevant problems. We publicly release our implementation through CaPTk (Cancer Imaging Phenomics Toolkit) and as an MITK plugin.

Fully automated detection of paramagnetic rims in multiple sclerosis lesions on 3T susceptibility-based MR imaging.

BACKGROUND AND PURPOSE: The presence of a paramagnetic rim around a white matter lesion has recently been shown to be a hallmark of a particular pathological type of multiple sclerosis lesion. Increased prevalence of these paramagnetic rim lesions is associated with a more severe disease course in MS, but manual identification is time-consuming. We present APRL, a method to automatically detect paramagnetic rim lesions on 3T T2*-phase images. METHODS: T1-weighted, T2-FLAIR, and T2*-phase MRI of the brain were collected at 3T for 20 subjects with MS. The images were then processed with automated lesion segmentation, lesion center detection, lesion labelling, and lesion-level radiomic feature extraction. A total of 951 lesions were identified, 113 (12%) of which contained a paramagnetic rim. We divided our data into a training set (16 patients, 753 lesions) and a testing set (4 patients, 198 lesions), fit a random forest classification model on the training set, and assessed our ability to classify paramagnetic rim lesions on the test set. RESULTS: The number of paramagnetic rim lesions per subject identified via our automated lesion labelling method was highly correlated with the gold standard count per subject, r = 0.86 (95% CI [0.68, 0.94]). The classification algorithm using radiomic features classified lesions with an area under the curve of 0.82 (95% CI [0.74, 0.92]). CONCLUSION: This study develops a fully automated technique, APRL, for the detection of paramagnetic rim lesions using standard T1 and FLAIR sequences and a T2*phase sequence obtained on 3T MR images.

Association between fine mapping thymic stromal lymphopoietin and atopic dermatitis onset and persistence.

BACKGROUND: Atopic dermatitis (AD) is a common chronic relapsing skin disease. Genetic variants have been associated with skin barrier function and immune regulation. Thymic stromal lymphopoietin (TSLP), an immune regulator, has been previously associated with AD. OBJECTIVE: To fine map TSLP and evaluate associations with the onset and persistence of AD. METHODS: TSLP variation was determined using targeted massively parallel sequencing in a longitudinal cohort of children with AD. Evaluations included linkage disequilibrium and the persistence of AD for as many as 10 years of follow-up. The association between the presence of AD and rs1898671 variation was evaluated in a second independent cohort. RESULTS: The minor variant frequency for rs1898671 was 23.5% (95% CI, 21.4%-25.8%). This variant was not in linkage disequilibrium with other TSLP variants in the longitudinal cohort (n = 741). White children with AD were less likely to have rs1898671 variant (odds ratio [OR], 1.41; 95% CI, 1.20-1.66) than Genome Aggregation Database controls. Children with AD and the rs1898671 variant during follow-up were more likely to have remission than children who were wild type for rs1898671 (OR, 1.56; 95% CI, 1.26-1.91). In the second cohort (n = 585), the rs1898671 variant was less prevalent in those with AD than those without. The protective effect was greater in rs1898671 heterozygotes (OR, 1.91; 95% CI, 1.34-2.75) than homozygotes (OR, 1.28; 95% CI, 0.61-2.70). CONCLUSION: TSLP and specifically rs1898671 are important in the pathogenesis of AD and could represent a potential clinical target for the development of therapies to treat individuals with AD.

Protein-structure-guided discovery of functional mutations across 19 cancer types.

Local concentrations of mutations are well known in human cancers. However, their three-dimensional spatial relationships in the encoded protein have yet to be systematically explored. We developed a computational tool, HotSpot3D, to identify such spatial hotspots (clusters) and to interpret the potential function of variants within them. We applied HotSpot3D to >4,400 TCGA tumors across 19 cancer types, discovering >6,000 intra- and intermolecular clusters, some of which showed tumor and/or tissue specificity. In addition, we identified 369 rare mutations in genes including TP53, PTEN, VHL, EGFR, and FBXW7 and 99 medium-recurrence mutations in genes such as RUNX1, MTOR, CA3, PI3, and PTPN11, all mapping within clusters having potential functional implications. As a proof of concept, we validated our predictions in EGFR using high-throughput phosphorylation data and cell-line-based experimental evaluation. Finally, mutation-drug cluster and network analysis predicted over 800 promising candidates for druggable mutations, raising new possibilities for designing personalized treatments for patients carrying specific mutations.